2D Rashba system in AC magnetic field 
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The response of an electron system to a DC measurement electric field has been investigated in the 
case when the system is driven out of the equilibrium by the magnetic ultra-high frequency field that 
leads to combined transitions. The discussed model includes contributions from Landau quantization 
and from microwave irradiation. Impurity centers are considered as sources of scattering. It has 
been shown that the perturbation of the electron system by the ultra-high frequency magnetic field 
leads to oscillations of the diagonal components of the conductivity tensor. 
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I. INTRODUCTION 

The interest in theoretical studies of transport phe- 
nomena in 2D electron systems has raised substantially 
after the discovery 0, of oscillations of the diagonal 
components of the conductivity tensor in "ultraclean" 
GaAs/Al^Gai-^As samples in the classical interval of 
the magnetic field intensity, where the Shubnikov — de 
Haas (SdH) oscillations don't manifest themselves. To- 
gether with the oscillations of the diagonal components of 
the conductivity tensor caused by absorption by charge 
carriers of the energy of the microwave radiation field 
and transitions between Landau levels, "beats" have been 
found experimentally in the interval of more weak 
magnetic fields. Such beats are usually related with 
the manifestations of the interaction between kinetic and 
spin degrees of freedom of the conductivity electrons. 
Such interaction is the spin-orbit interaction (SOI) that 
is known to be the origin of numerous effects in transport 
phenomena observed in such systems. Among them there 
are, e.g., beats in SdH oscillations spin accumulation 
, magneto-electric effect , etc. SOI also leads to the 
possibility of electron transitions between Landau levels 
in the magnetic field at the combined resonance frequen- 
cies @ , thus transitions being possible both in antinodes 
of electric and magnetic fields Q ■ Finally, operation of 
a spin transistor (schemes of which have been considered 
in |3) is based upon spin degrees of freedom. All of the 
above has determined elevated interest to investigations 
of the SOI in 2D semiconductor structures. 

For the purpose of studying the SOI, it appears in- 
teresting to investigate a model in which the role of the 
SOI should manifest itself the strongest. Since the SOI 
depends upon both translational and spin degrees of free- 
dom, then it is a channel over which energy (both elec- 
tric and magnetic) can be absorbed from the ultra-high 
frequency field, thus causing transitions between Lan- 
dau levels. Because of that, it is interesting to investi- 
gate the response of a non-equilibrium electron system 
to a DC weak ( "measurement" ) electric field for the case 
when the initial non-equilibrium state is created with a 
high-frequency AC magnetic field that leads to combined 



transitions. The question is how this perturbation affects 
transport coefficients, in particular, the conductivity ten- 
sor. 

The discussed model includes the contributions fom 
Landau quantization and (in th long-wavelength limit) 
from the microwave radiation exactly, without use of the 
perturbation theory. We consider impurity centers for 
the role of scatterers, treating the scattering process per- 
turbatively. 



II. EFFECTIVE HAMILTONIAN 

The Hamiltonian of the system under consideration is: 

H(t) =H k + H s + H ks + H eh (t) + H° ef +H V + H ev . (1) 

Here H — k and H s are kinetic and Zeeman energies, 
respectively, in the magnetic field H = (0, 0, H): 



(Pi - {e/c)A{ Xl )f 
2m 

H s = hw s Y, Sf, hu s = gf i H. (2) 



Sf and pf are operators of the components of the 
spin and kinetic momentum of the ith electron, where 
\p a iPj] — —SijimhLu c s a /3 z , lu c — \e\H/mc is the cyclotron 
frequency, and /io is Bohr magneton. Hjjf is the hamil- 
tontal of the electrons' interaction with the electric field 
E = (E X , 0,0): 



H, 



ef 



(3) 



H e h(t) is the interaction of electrons with the AC mag- 
netic field: 



H eh (t)=gfi H(t)^2Si 



(4) 



H{t) = (H x (t),H y (t),H z (t)). Hev and H v are Hamil- 
tonians of the electron-lattice interaction and of the lat- 
tice itself, respectively. H ks (j>) is the interaction between 
translational and spin degrees of freedom. Its most gen- 
eral form is: 
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H ks (p) 



(5) 



2 



Here f(pj) is a pseudo- vector, components of which can 
be represented as a form of order s in the components of 
the kinetic momentum pj. 

The spin-orbit interaction leads to correlation of spa- 
tial and spin motion of electrons, thus, the translational 
and spin-related subsystems are not well-defined. Since 
the SOI is in some sense small, then one can perform a 
momentum-dependent canonical transformation that de- 
couples kinetic and spin degrees of freedom. Naturally, 
all other terms in the Hamiltonian, describing the inter- 
action of electrons with the lattice and external fields (if 
any) also undergo the transformation. In this case, the 
effective interaction of electrons in the system with exter- 
nal fields appears, which leads to resonant absorption of 
the field energy not only at the frequency of the paramag- 
netic resonance uj s or cyclotron resonance u> c , but also at 
their linear combinations, i.e. the combined resonance. 
The gauge-invariant theory of the combined resonance 
has been developed in |lCj . 

Assuming the SOI to be small, we perform the canon- 
ical transformation of the Hamiltonian. Up to the terms 
linear in T(t), we have: 



H = e T{p) He~ T{p) « H+ [T(p),H]. 



(6) 



The operator of the canonical transformation T(p) has 
to be determined from the requirement that, after the 
transformation, the k and s subsystems become indepen- 
dent. This requirement can be written as the following 
condition: 



Obviously, the criteria for the applicability of this theory 
is that, for characteristic values of the electron momen- 
tum p, the inequality ap <C h(u> c — u> s ) should hold. 

One can write the transformed Hamiltonian in the fol- 
lowing form: 

H(t) = H + H° ef + H eh (t) + [T{p), H eh (t) + 

+ H Q ef +H ev ]. (10) 

Hq = H k + H s + H v + H ev . 

Using the explicit expression for the operator T(t), we 
find: 



gf x [T(p),S a }H a (t) = 



2h(u; c - ui s ) 
T z+ H-(t)) + (T- + -T+-)H z {t)}. (11) 



Now we find the operators of power H i n l \if) = 
{ih)- l [Hi,H eh {t) + [T(p),H eh (t)]\ absorbed by kinetic 
(i = k) and spin (i — s) subsystems due to the inter- 
action of electrons with the AC magnetic field. We have: 



H k (h){t) 



gap Q uj c 



{(T~+ +T+-)H z {t)- 



H ks {p) + [T(p),H k +H s ]=0. 



(7) 



2h(u> c - LU S ) 

-(T z -H+(t)+T z+ H-{t))}, (12) 



Note that, after the canonical transformation, the H k 
and H s operators are the integrals of motion if H ev = 
and there is no interaction with external fields. 

We assume the specific form of the SOI term, namely, 
Rashba interaction, which is non-zero even in the linear 
order in momentum: 



H ks (jp) = ae ztk ]T S[,/: = ^ ^(Sfpj - S7p+), (8) 



Hs(h)(t) 



S-H+(t)-S+H-(t))- 



-(T-+ +T+-)H z (t) (13) 



2h(w c - w a ) 
The total absorbed power can be written as: 

H k(h) (t)+H s{h) (t) = jgHP(t), (14) 



S ± = S x ±iS\ p ± =p x ±ip v . 

Here a is the constant characterizing the SOI, e is the 
fully-antisymmetric Levi — Chivita tensor. 

Now we find the explicit expression for the operator 
T(p). Inserting the operator (JSJ into the general solution 
for Eq. JjJ and integrating, we obtain: 



u 

T(p) = ~ lim [ dte £t e ltH ° /h H ks ip)e- UH ° /h 

h e^+o j 



< J >>) 



2h(u; - u s ) 



J2^pj- s Jpi)- ( 9 ) 



9^0 
ih 



\S^[Tip),S%H k ip)+H s 



The interaction of the spin degrees of freedom of the 
conductivity electrons with the AC magnetic field H e h [t) 
leads to resonant transitions at the frequency oj s . How- 
ever, as one can see from the expressions above, the effec- 
tive interaction [Tip), H e h(t)] leads to combined transi- 
tions at frequencies lo c ± u> s and the cyclotron frequency 
lu c . Since, for our further calculations, the response of 
the non-equilibrium system to the measurement electric 
field is interesting, in which the contribution from the 
translational degrees of freedom dominates, we will re- 
strict our consideration to the effective interaction solely. 
Besides that, we limit the consideration to the case when 
the DC and AC magnetic fields are perpendicular to each 
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other: H(t) = (H x (t), H y (t),0). In this case, the effec- 
tive interaction responsible for the combined transitions 
has the following form: 

H eh ,i{t) = [Tip), H eh (t)] = 



2H(uj c - uj s ) 



H SZ M 



+ —iuit 



p-e^). (15) 



lo\ s = geHi/ (2mQc). Hi is the intensity of the circularly 
polarized magnetic field, rotating with the frequency to. 

The dependence of the effective interaction H e h,i(t) 
upon time causes certain difficulties while calculating 
the non-equilibrium response of the electron system to 
the measurement electric field. Thus, it is expedient to 
carry out one more canonical transformation (Appendix 
A), that removes the interaction H e h,\{t) and renormal- 
izes the electron-impurity interaction Hamiltonian (Ap- 
pendix B), which acquires the time dependence then. In 
the canonically transformed system, impurities act as a 
coherent oscillating field that leads to resonant transi- 
tions. 



III. NON-EQUILIBRIUM RESPONSE 

We assume that the initial non-equilibrium state of the 
system under consideration is created by the ultra-high 
frequency magnetic field and can be described with the 
distribution p(t). Obviously, if some additional pertur- 
bation acts upon the system, then a new non-equilibrium 
state is formed in the system, that requires an extended 
set of basis operators for its description. The new non- 
equilibrium distribution is described with the operator 
p(t,0). The task is to find the response of a non- 
equilibrium system to a weak measurement field. 

We write the operator p(t) using the integral represen- 
tation for the non-equilibrium statistical operator 8] . In 
the linear approximation in the external field E, we have: 



p(t)=p(t)-i J dt l e etl U{t + t l )x 

— OO 

xiL° efP (t + ti)U + {t + tx). (16) 

iL° ef A=^[A,H° ef ]. 
The operator p(t) satisfies the equation 

^ + ±{p(t),H(t)} = -e(p(t) - Pq (t))> (17) 
which has the following solution: 

o 

p(t) = p q (t)- / dt ie etl u(t + h)x 



iL(t)A=-[A,H(t)}, 
in 

where U(ti) is the evolution operator that is in fact a 
chronologically ordeted exponent and satisfies the differ- 
ential equation 



dU(t) 
dt 



in 



The quasiequilibrium statistical operator p q (t) can be ex- 
pressed as: 



pS) 



-S(t) 



S(t)=$(t)+Y,P+F n (t), (19) 



where S(t) is the entropy operator, $ is the Massieu — 
Planck functional. P n , F n are sets of the basis operators 
and their conjugate functions, that describe the electron 
system. 

Characterizing the state of the system by the mean 
values of the operators p, H s , N, H v (N being the 
operator of the number of electrons), we obtain for the 
entropy operator: 

S(t) = + p k (H k - V(t)p- p'N) + (5 S H S + 

+ (3H V = S° + AS(t). (20) 

AS(t) = -f3 k V(t)p. 

Here Pk, Ps, p' — P — mV 2 /2, V, P are parameters ther- 
modynamically conjugate to the average values of the in- 
troduced operators. They have the meaning of inverse ef- 
fective temperatures of the kinetic and spin electron sub- 
systems, chemical potential, drift velocity and the inverse 
lattice temperature. Introduction of the effective temper- 
atures allows one to treat effects related to "heating" of 
the electron and spin subsystems by external fields. 

A linear addmittance corresponding to an arbitrary 
operator B in the case when the external force oscillates 
harmonically in time with the frequency lo\ can be ex- 
pressed as: 



XBA(t,UJl) 



J in 



x^ L [A,p{t + t x M} (21) 



Within the framework described above, the task of 
the admittance calculation is reduced to obtaining the 
transport matrix Tba(^^i), which plays in the non- 
equilibrium case the same role as in the case of the equi- 
librium response: 

XBA{t,ui) =Xs4t,0)- — 77 r— — ■ (22) 

T B A{t,uJi) + e - iuji 



x {— +iL(t)} Pq (t + t 1 )U+{t + t 1 ), (18) 



XBA(t,0) = (B,A), 



T, 



BA 



(B,A)» 



-{B,A) W1 . (23) 
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Here 



u 

(B,A) = ~ J dhe^ Sp{Be ltli [Ap(i + ii,0)]}, 



(24) 



(-B.A)"' 1 =-i [ dt 1 e {£ -^ l)tl dt 2 e et2 Sp{Bx 
in J 

— CO 

Xe i < tl+ta > x Kp(t + ti+t 3 ,0)]}. (25) 

The real part of the transport matrix determines the 
relaxation frequency of the non-equilibrium electrons' 
momentum. 



IV. MOMENTUM RELAXATION RATE 

Assuming the temperatures of the translational and 
spin subsystems to be equal (that corresponds to neglect- 
ing any heating effects), in the Born approximation upon 
the electron-scatterer interaction we obtain for the relax- 
ation frequency: 

n l 

1 (3 1 



t 2mn 



Re ih J *i e(£_Wl)tl *2e £t2 J dX > 

— oo 



+ t 2 ),H k +H s y q -*}, p^ = -\p a ,H ev ]. (26) 

Now we expand the formula (|26[l using the explicit ex- 
pression for the renormalized electron-impurity interac- 
tion H ev from Appendix B. 

It is convenient to introduce the following notation: 

A(\,t 1 + t 2 ) = sp{p+ ) (ty L <>^ +t ^x 

x P x q [p^(t + t 1 +t 2 ),H k + H s ] Pq - x . (27) 

Inserting the explicit expression for the electron-impurity 
interaction and averaging over the system of scatterers, 
we obtain: 

A(X,t 1+ t 2 ) = ]T \v(q)\ 2 N ie ~ u ^ + ^jH\K q \)q 2 x 

in' i 

x Sp{25;e^e 4io( * 1+t2 V?[25*|,e-^' , H k +H s ]p\- X }. 

(28) 

Here Ni is the impurity concentration. Rewriting this 
expression using secondary quantization and averaging 
the electron operators using Wick's theorem, we have: 

A(X,h+t 2 ) = iViq^N.e-^+^x 

x e -i( t i+ t 2)(£M-<v) e "' 3 <=( e f -<v) A x 
x jf{\K q \) q \e v - £ai )|(2SV^W| 2 /,(1 - U), (29) 



where / is the Fermi — Dirac distribution. 

After integration over A, t\ and t 2 , this produces: 



dhe {s - luJl)tl dt 2 e Et2 / dXA^h+h) 



= J2 |^(g)| 2 ^Jf(|^|)< Z 2 |(25 2 e^)^| 2 (/, - /, 

e — iuii — Uuj — (i/h)(e^ — e v ) 
In the e — > limit, we obtain: 



(30) 



m- 1 



e — ilu)— {i/K){£^ - £ v ) 
1 



mSQfiw + ( e -e w ))), (31) 



Ihw + {Ep - £„) 

where P denotes the Cauchy principal value. Taking the 
u>i — > limit (because calculating a zero-frequency re- 
sponse is our goal), we have: 

qfivl 

x \(2S z e^ r U\ 2 (f(£ + lhw)- f{£))x 
d_ 

We 

The equation 1)32(1 contains a singularity in its right 
hand side, which is removed, as usual, due to broadening 
of the Landau levels by scattering electrons on impurities: 



x S(e-£^)— 5(lhu + E -£ v ) (32) 



^72 



• exp 



2r 2 



(33) 



The Landau level width T can be expressed via the elec- 
tron mobility a in zero magnetic field: 



T = h 



2j n uj c mu 

, T tr = -rj (34) 

ITTtr e 



Note that, if T > T, one can pull f(£ ± hto) - f(£) out 
of the £ integral as a slowly-changing factor. This yields: 



J d£-^D v {£±hw)D^{£) 



4r 3 



exp 



ijj, - £ v ± hmf 

4r 2 



(35) 



The wave functions upon which the matrix elements in 
(|32|l are calculated, have the form: 



i> v = i>nk*s* 



1 



x exp 



V2 n n\ir 1 / 2 £ 

(y - yo) 2 



2P 



exp(ik x x)x 

H n (^-) X s* (36) 
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yo = i 2 k x is the cyclotron orbit center coordinate, £ is 
the magnetic length, H n (x) denotes Hermite polynomi- 
als, and xs z is the eigenfunction of the z spin projection. 
Calculation of the matrix element in 1|32[1 produces: 



\(v\2S z exp(iqr)\ti)\ 2 = cxp 

(minfj^, ri^))! / £ 2 q 



fq 2 



(max(n„,n M ))! 



L 



\nv—n»\ 



S SS,S*Jk*,q x +k* (37) 



In the case of sufficiently weak AC magnetic field one 
can neglect terms with |Z| > 1 and use the approximation 
J±i(x) = ±x/2. As a result of such simplifications, we 
obtain: 



1 _ h 

V ~~ 8mn£ 2 



IViq^NtlK^q 2 

qn v n^l— ±1 



x exp 



^ 2 i7 2 \ (min(n I/ ,n M ))! /^ 2 g 



(max(n I/ ,n Jli ))! 



2„2 



(^-)) (/M-/(e M ))x 



3 / 2 (e ; . - 



ihw) 



4r 3 



■ exp 



(e M — e v + Zftw) 1 



4r 2 



(38) 



Integrating over q, we have: 



d(q 2 )q exp 



t 2 q 2 
l 2 q 2 



l 2 q 2 



2 

)) = 



8 (max(n„,)i M ))! 
£ 6 (min(n,/,7fy))! 



x (n 2 + n 2 + 3(71,, + n M ) + ^ + 2) (39) 

Thus, for the case of point scatterers, where V(q) does 
not depend on q, the radiation-induced correction to the 
inverse relaxation time is: 



1 _ h 

t Amni 8 



E 



(lO c - L} s ) 2 (uj - LO c ) 2 



x (n 2 + n 2 + 3(n„ + n M ) + 4n„n jU + 2) x 



X (/(£„) -/(£„))- 



rV2 



ihw) 



x exp 



r 3 

(e M - £y + ^w) 2 
if 2 



(40) 



Using the expression for the momentum relaxation rate, 
one can also write the formula for the diagonal com- 
ponetns of the conductivity tensor <j xx , according to 
which the numerical calculations in the next section are 
carried out: 



(41) 



m u> 2 + t 
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FIG. 1: Photoconductivity of the 2D electron gas vs of the 
magnetic field induction for different values of electron mo- 
bility. The radiation frequency is 50 GHz and 7 = 2. 



V. NUMERICAL ANALYSIS 

Numerical calculations of the diagonal components of 
the conductivity tensor according to Eq. I|41|l have been 
carried out with the following parameters: m = 0.067m 
(mo is the free electron mass), the Fermi energy is £p = 
10 meV, the mobility of the 2D electrons varies as \i f=s 
0.1 — 1.0 x 10 7 cm 2 /Vs, the electron density n = 3 x 
10 11 cm~ 2 . The microwave radiation frequency is / = 
50 GHz, the temperature is T w 2.4 K. The magnetic 
field varied as 0.02 - 0.3 T. 

The dependence of the 2D electron gas photoconduc- 
tivity on the oj/uJc ration is presented in Fig.^ One can 
see that the dependence of electron mobility upon the 
magnetic field has the oscillating character. In the region 
of low magnetic field, the oscilation amplitude drops sig- 
nificantly when zero-magnetic-ficld mobility is decreased. 

In Fig. the photoconductivity dependence upon the 
magnetic field is presented for different values of 7„ and 
the same microwave radiation frequency 50 GHz and elec- 
tron mobility /i = 0.6 x 10 7 cm 2 /Vs. As one can see, 
the oscillation amplitude is very sensitive to the width of 
Landau levels. 



VI. CONCLUSION 

The response of a non-equilibrium electron system to 
the DC electric measurement field has been studied for 
the case when the initial non-equilibrium state of the 
system is created by an ultra-high frequency magnetic 
field that leads to combined transitions. Within the pro- 
posed theory, it has been shown that such perturbation 
of the electron system essentially influences the trans- 
port coefficients and leads to the oscillations of the di- 
agonal components of the conductivity tensor. The dis- 
cussed effect is analogous to the phenomenon observed 
in GaAs/AlGaAs heterostructures with ultra-high elec- 
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FIG. 2: Photoconductivity of the 2D electron gas vs of the 
magnetic field induction for different values of electron mo- 
bility. The radiation frequency is 50 GHz and 7 = 2. 



tron mobility [2j. However, unlike that phenomenon, the 
manifestation of the oscillatory pattern is dictated by the 
spin-orbit interaction existing in the crystals under con- 
sideration. 



APPENDIX A 

In this Appendix, a canonical transformation W2(t) is 
built, that excludes the renormalized interaction with the 
AC magnetic field from the effective Hamiltonian. The 
canonical transformation operator is searched from the 
equation: 



tW^(t)(~ih^- t +H k + H s + H ehA {t))W 2 (t) = 

= -ih%- + H k + H a . (Al) 
ot 

The operator W 2 if) is searched in the following form: 



W 2 (t) = exp(*£>-(«)p+^ + v + (t) P jS*))x 
3 

x exp(i0(i)), (A2) 

where one has to determine the parameters 9(t) and 

In order to determine these parameters, the canonical 
transformation is applied to all terms in the left-hand 
side of Eq. (33): 



8 



8 



Wl{t)(-ih—)W 2 (t) = -ih— + M{t)+ 



dt 



Ot 



Kj2(v-(t)pt+V + (t)pj) (A3) 



4imhuj c T] + (t)SjPj + ^imhhj c rj (i)SjPj 

+ m 2 h 2 uj 2 c ri + (t)r 1 -(t)) (A4) 



W 2 \t)H s W 2 {t) = H s 



(A5) 



W%(t)H eh>1 (t)W 2 (t) = 



iauji s 



2(lu c - LU S ) 



- 2imhw c r 1 + {t)S*)e- iuJt - 

- (pj + 2imhw c rf (t)S?)e iut ) (A6) 

Obviously, the equation (|A1|) becomes an identity un- 
der the following conditions: 



lawis ±i 



2(uj c - u s ) 



0, (A7) 



. • , . Nmh 2 uj 2 , , . , . 

m) + — r ^v + (t)v-(t)+ 

+ -r r Jf r * e + r,-(ty^) = 0. (A8) 

Thus one can write down the explicit expressions for 



6(t) 



aujis 



2H(uj c - lo s )(lu - LO c ) 



Nma 2 Lol s {2,uj c - Au)t 
16H(u) c — lu s ) 2 (lu — u> c ) 2 



(A9) 



(A10) 



Therefore, the explicit form of the canonical transfor- 
mation W2 is known. The physical meaning of this trans- 
formation is the change to two non-uniformly translation- 
ally moving reference frames, different for electrons with 
opposite spin directions. 



APPENDIX B 

Since we consider only elastic scattering, the electron- 
impurity interaction Hamiltonian has the form: 



(Bl) 



where V(q) is a Fourier component of the potential cre- 
ated by a single impurity corresponding to the wave vec- 
tor q, p(q) is a Fourier component of the impurity num- 
ber density. 
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As a result of the canonical transformation \V2if), _ auJ is{Qx — i%) 
renormalization of the electron- impurity interaction hap- q {u> c — u s )(ui — ui c ) 
pens. For obtaining the renormalized Hamiltonian of the 
electron-impurity interaction, it is sufficient to calculate 
W\ (t) exp(iqrj)W2(t). Using the explicit form of the op- 
erator W|(t), we obtain: Thus, the renormalized Hamiltonian of the electron- 



w\(t)e iqr ^W 2 (t) = exp{iq( rj + Arj)), (B2) 

where 



impurity interaction is: 



Ar4 = —, ^ -S? x (coswt,sinart,0) (B3) ^ 



functions J t (x) yields: x (2S'--^ T e^M Ji{\K„\). (B5) 

* i\KJ 1 



Expansion of the exponent in (|_B2|I in terms of Bessel 
nctions Ji{x) yields: 

e -iR.(2jc,s Je *") = g f 2 s;-^-e^\ miks, 

I — — 00 ^ I Ql / 

(B4) 
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